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| Abstract. We report results from Exploratory Time observations of 

\£> • the double-pulsar system PSR J0737— 3039 using the Green Bank Tele- 

scope (GBT). The large gain of the GBT, the diversity of the pulsar 
backends, and the four different frequency bands used have allowed us to 
make interesting measurements of a wide variety of phenomena. Here we 
briefly describe results from high-precision timing, polarization, eclipse, 
scintillation velocity, and single-pulse work. 
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, 1. Introduction 

; 

In December 2003, we proposed for and were awarded Exploratory Time to 
observe the spectacular double-pulsar system PSR J0737— 3039 (hereafter 0737; 
+3 , Burgay et al., 2003; Lyne et al., 2004, and contributions to this volume from 

Burgay, Manchester, Kramer, and others) as part of the NRAO Rapid Response 
Science program 1 . We observed 0737 five times at four different frequencies 
(427, 2x820, 1400, and 2200 MHz) using the Berkeley-Caltech Pulsar Machine 
(BCPM) and the Green Bank Pulsar Processor (GBPP; which measures full 
polarization information). For two of the observations (one each at 427 and 
820 MHz) we also used the new GBT Spectrometer SPIGOT card (a correlator- 
based instrument that outputs lags at 25MB/s) in some of its first scientific 
observations. These data are public and can be obtained from NRAO. 



2. High-precision Timing 

The high sensitivity of the GBT allowed us to measure high signal-to-noise 
ratio (SNR) pulse profiles from PSR J0737-3039A (hereafter A, see Fig. 1) 
at each observing frequency, precise (^11 /is) pulse arrival times (TOAs) from 
A in 1— 2min integrations, and TOAs throughout ~80% of the orbit for PSR 
J0737— 3039B (hereafter B). Without a long timing baseline we could only mea- 
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sure "local" parameters, such as the Shapiro delay "shape" s = sinz for A, and 
better projected orbital semi-major axes x = asini/c for the pulsars. 

For A we used the Lyne et al. (2004) ephemeris as our base and extracted 
TOAs from the high SNR 820- and 1400-MHz BCPM and SPIGOT observations. 
Using TEMPO 2 , we determined s=0.99962+o o[Jq95 (corresponding to z=88?4±^4) 
and ^=1.4150342(82) lt-sec with a post-fit RMS of ~11 fis (the SPIGOT data 
alone had an RMS of <^7/Us). The errors for s and xa were estimated using a 
bootstrap analysis. We note that we had to fit out apparent linear arrival time 
drifts of order tens of /us during each of the observations. We believe that these 
(likely) systematic drifts could result from uncorrected variations in the linearly 
polarized portion of the pulse profile due to the rotation of the receiver feed with 
respect to the sky and gain differences between the two orthogonally polarized 
IFs. See Ransom et al. (2004) for more details of the A timing. 

To time B, we determined a set of 29 template profiles (all referenced to 
the same gaussian component) as a function of the orbital phase (since B's 
pulse profile changes systematically throughout the orbit) and extracted TOAs 
from each of the BCPM observations. Using the Lyne et al. (2004) timing 
solution as a starting point, we used TEMPO to fit only for xb and the local spin 
frequency. The post-fit residuals show significant systematics as a function of 
orbital phase (see Fig. 3), implying that there may be no stable timing reference 
point in the varying profile (i.e. A's wind likely affects the radio beam of B 
as a function of orbital phase). If we multiply TEMPO'S formal errors on xb by 
10.3 (the square-root of the reduced-^ 2 of the fit) to try and account for some of 
these systematics, we get xb=1. 5126(14) lt-sec and a mass ratio R = mA/mB = 
xb/xa = 1.06893(97), which is a factor of ~6 more precise than reported in 
Lyne et al. (2004). A Ma vs. Mb diagram based on the GBT data but using 
the reported value of u> from Lyne et al. (2004) is shown in Fig. 1. 

3. Polarization 

The GBPP coherently dedisperses and folds full-Stokes data on a single pulsar in 
real time. The most useful polarimetric data sets were those taken at 820 MHz: 
B on 2003 Dec 19 and A on 2003 Dec 24. We used the full track on A to perform 
polarization self calibration of the data. This calibration fit included terms 
accounting for channel gain difference, feed cross coupling, feed ellipticity, and 
cable delays. The RMS average of all such off-diagonal Mueller matrix elements 
was 0.085, so we are confident that our calibration is fairly good. We also used 
measurements of the observatory correlated noise source and observations of 
PSR B1929+10 to confirm our calibration. 

Fig. 4 shows the 820 MHz GBPP observations of A and the mirror-like 
symmetry of the pulse profile. We have chosen longitude (ft = as the center of 
this symmetry, which is visible in both the total power and linearly polarized 
profiles. Also shown is the position angle versus pulse longitude for points having 
>4o" significance. The longitude region —100° < (ft < —70° appears to contain 
an orthogonal emission mode jump. The region +90° < (ft < +105° may also be 
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Figure 1. (Left) Average pulse profiles of A from the GBT. Changes 
in the relative pulse amplitudes with observing frequency are obvious. 
The short lines above and to the right of the profiles depict the time 
resolution of the data. (Right) A diagram of the Ma vs. Mb plane 
(after Lyne et al. , 2004) based on the GBT data (except for the advance 
of periastron Cj). The shaded region requires an inclination i > 90° 
and is therefore not allowed. The other plotted parameters are the 
mass ratio R = m/i/mg, the Shapiro delay "shape" parameter s = 
sini, and the i as determined using both timing (§2.) and scintillation 
measurements (§5.). Lower and upper limits on the various parameters 
are labeled with LL or UL respectively. 
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Figure 2. Timing residuals from the 820 MHz SPIGOT observation 
of A before (top) and after (bottom) the addition of a Shapiro delay 
"shape" fit for a 1.25 Mq companion. The Shapiro delay is easily 
detected during each orbit. The post-fit RMS residuals for these data 
alone are ^7/is. 
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Figure 3. Post-fit timing residuals as a function of orbital phase for 
B. The only fit parameters were the projected semi-major axis xb and 
the spin frequency. The TOAs are from BCPM data from all five 
GBT observations of 0737. Systematic drifts in the TOAs are obvious, 
even after using orbital-phase-dependent profile templates to extract 
the TOAs. The "bright" portions of the orbit as defined in Demorest 
et al. (2004) and Ramachandran et al. (2004) are labeled as I, II, III, 
and IV. The possible "eclipse" of B is also indicated (see §6.). 



affected by quasi-orthogonal modes. We find the rotation measure of this pulsar 
to be -112.3 ± 1.5radm~ 2 . 

The mirror symmetry of the profile suggests that the emission is coming 
from two traversals of a wide, hollow emission cone, rather than from two op- 
posite magnetic poles. We have performed a rotating vector model (RVM; Rad- 
hakrishnan & Cooke, 1969) fit to the data which seems to support this interpre- 
tation. The fit gives values for the angle between the magnetic and rotation axes 
(a) and the closest angle between the magnetic axis and the line of sight ([3). We 
can also calculate the cone opening half-angle p, using the longitude separation 
of the two peaks. The fit produces two solutions, one with a = 4° ± 3° and one 
with a = 90° ± 10°. The large-a solution requires a extremely large opening 
angle, p ~ 90°. We favor the small-a solution, as it is more consistent with the 
mirror symmetry, and permits more modest values of p. In this solution, (3 is 
unconstrained, although we argue that it is likely that (3 ~ p. 

The polarimetric properties of B are also interesting, though less easily 
understood. In addition to showing a remarkable intensity variation through- 
out its orbit, this pulsar shows variation in its degree of linear polarization. 
The brightest "window" (I), just before conjunction, shows a peak linear polar- 
ization fraction of ~14%. However, in the next bright window (II), just after 
conjunction, the pulse is almost completely unpolarized (^1%). From window 
I, we determine the rotation measure to be —118 ± 12radm~ 2 , consistent with 
our value for A. More details on these observations, the model fits, and their 
implications are given in Demorest et al. (2004). 
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Figure 4. (Top) GBPP 820 MHz profiles of A. The dotted and 
dashed lines show linear and circular polarization respectively, while 
the solid line is the total intensity. The mirror-like symmetry of the 
profile (see inset) implies that the pulses come from two crossings of a 
single cone of emission. (Bottom) A rotating vector model (RVM) fit 
to the data (see §3.). See Demorest et al. (2004) for more details. 



4. Eclipses of A 

We used the BCPM data to study the eclipse of A. We folded data obtained 
within 2min of A's superior conjunction in 2-s intervals using 256 phase bins 
and then cross-correlated each profile with a template. The cross-correlation, 
done in the Fourier domain, assumes the data profile p has the form p(j) = 
a + b x s(j — t) + n(j), for 1 < j < 256, and where s(j) is the template, b is a 
scale factor, n(j) is the noise background, and r is a temporal offset. Values of 
a, b and r were adjusted iteratively to minimize a x 2 fit to each profile (Taylor, 
1990). We used the scale factor b and its standard deviation as a relative measure 
of the pulsar's flux and its uncertainty, respectively. Our resulting eclipse light 
curves at three radio frequencies are presented in Fig. 5. The eclipse properties 
are remarkably independent of frequency, as previously noted qualitatively by 
Lyne et al. (2004). Because of its largely achromatic properties, we show in 
Fig. 5d an eclipse light curve obtained by averaging those in panels a, b, and c. 
A striking eclipse property is that it is asymmetric at all frequencies in a similar 
way. The eclipse ingress is clearly longer in duration than the egress. 

To quantify some of these properties, we fit the light curves with two func- 
tions of the form F(<j)) = (e^-^o)/™ + ij-i^ w here F(<j>) is the flux at orbital 
phase (p, defined such that superior conjunction is at = 0. The best fits are 
shown as solid curves in Fig. 5 and all yield values of reduced-x 2 near unity. 
Using the fit parameters, we quantitatively show that the eclipse duration is 
independent of frequency with a limit on a power-law index of ~0.5; that the 
ingress takes ~3 times longer than egress, and that the eclipse very likely lasts 
longer post-conjunction. There is a hint that the eclipse asymmetry grows with 
increasing radio frequency, but this is not statistically significant. 
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Figure 5. Plots of the eclipses of A at 427-, 820- and 1400-MHz with 
a time resolution of ~2s. The asymmetric ingress and egress of the 
eclipses and its near frequency independence are obvious. The deepest 
(i.e. most complete) part of the eclipse occurs just after conjunction. 



A possible model for the A eclipse is presented in these proceedings by 
Arons et al. (see also Arons et al. in prep). A's wind likely confines B's magne- 
tosphere on the side facing A, causing it to resemble a time-dependent variant 
of the Earth's magnetosphere. Synchrotron absorption in a shock heated "mag- 
netosheath" surrounding and containing B's magnetosphere is then responsible 
for the eclipse. Additional details on the data analysis and the eclipse model 
can be found in Kaspi et al. (2004). 



5. Scintillation Velocity 

The high SNR of the GBT data at 1400 MHz and the fantastic frequency res- 
olution of the SPIGOT at 820 MHz allowed us to make "self-calibrating" mea- 
surements of the orbital variation of the diffractive scintillation timescale Atj 
for A and hence, the so-called scintillation velocity Viss oc 1/Atd (e.g. Cordes 
& Rickett, 1998). Closely following the methodology of Ord, Bailes, &: van 
Straten (2002), we computed dynamic spectra for A for the two observations 
(see Fig. 6), and then using the known (from timing) orbital motion of A, fit a 
model consisting of four parameters: the orbital inclination i, the systemic veloc- 
ity parallel- (V^ane) and perpendicular-to (V perp ) the line of nodes on the sky, and 
the scaling parameter k. Our results indicate that 0737 is moving across the sky 
at Viss— 140.9±6.2 kms -1 , with component veocities of V^i ane ^96.0±3.7kms -1 
and V per p~103.1±7.7kms~ 1 . When combined with the timing fits described in 
§2., the best-fit inclination is i = 88?7 ± 0?9. 
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Figure 6. Dynamic spectra (top) and a scintillation velocity model 
fit (bottom) for the 820-MHz SPIGOT observation. Only the top 1/4 
of the band is plotted. The fit implies that 0737 is moving across the 
sky at ~140kms~ 1 . See §5. and Ransom et al. (2004) for more details. 



Using a simple kinematic argument and estimates of the characteristics of 
0737 before the second supernova (Dewi & van den Heuvel, 2004; Willems & 
Kalogera, 2004), the large V perp implies that B was born with a kick speed of 
^100 km s^ 1 . Future measurements of the scintillation velocity should allow us 
to remove the contaminating effects of the Earth's motion around the Sun and 
the differential rotation of the Galaxy. A precise systemic velocity in combina- 
tion with a VLBA-determined proper motion will provide an accurate geometric 
distance to the system (which could be crucial for high precision tests of gravi- 
tational theories; see the contribution by Kramer in this volume). See Ransom 
et al. (2004) for more details. 



6. Pulse Fluctuations and Morphology 

Due to the superior sensitivity of the GBT, we detected B at good signal-to-noise 
levels almost throughout the orbit. Fig. 7 shows the normalized flux received 
from B as a function of the true anomaly (referenced from the line of nodes) 
or b, at 820- (filled circles) and 1400-MHz (open circles). The light curve shows 
significant systematic variation, as first shown by Lyne et al. (2004), and four 
prominent phase windows where the pulsar is bright (indicated by I, II, III & IV; 
see also Fig. 3). In the range of 4> or b ~ 6° to 65°, the pulsar appears to exhibit 
what may be interpreted as an eclipse, where we did not detect any signal within 
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Figure 7. 820- (filled circles) and 1400-MHz (open circles) orbital 
lightcurve of B. See §6. and Ramachandran et al. (2004) for details. 



our sensitivity limits. In fact, around the eclipse boundary, the reduction in SNR 
was so significant that the pulsar became undetectable within a few rotation 
periods of B. The Arons et al. magnetopause model (these proceedings; 2004, 
in prep) may be able to explain at least some of B's systematic flux variation 
(including the eclipse), as well as A's eclipse. Jenet & Ransom (2004, see also the 
contribution from Jenet in this volume) propose an alternate model for B's flux 
variation where B brightens when it is illuminated by the cone of (presumably 
high-energy or particle) emission from A. 

Both A and B exhibit significant intrinsic pulse to pulse variations. The 
measured modulation index m at 820 MHz as a function of pulse longitude for 
A is in the range 0.6—2.0. B's profile structure is much more complex, where 
even the average profile changes systematically as a function of or b (as first 
discussed by Lyne et al., 2004), presumably indicating a change in the emission 
beam morphology. In windows I and II, values of m are in the range 1—1.5 and 
1—2, respectively. These values are significantly greater than what is expected 
from interstellar scattering (~7%), given the measured frequency and temporal 
dimensions of scintles and their filling factor (Ransom et al., 2004). The spectral 
shape of these intrinsic fluctuations are "noiselike" (i.e. white), indicating the 
lack of any organized drifting patterns. We present a much more detailed study 
of the above mentioned properties in Ramachandran et al. (2004). 



7. Conclusion 

The fantastic and highly-varied science to come from these five observations 
gives some indication of how important a role the GBT will play in future work 
on 0737. These observations provide the best constraints yet available on the 
inclination of the orbits, the mass ratio of the pulsars, the geometry of A's radio 
emission, the asymmetries and achromaticity of A's eclipses (and the nature of 
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B's magnetosphere that cause them), the orbital modulation of B's pulsed flux, 
and the systemic velocity of 0737 and the kick that caused it. 

The high gain of the telescope, the availability of many sensitive and useful 
receivers, and the new pulsar backends (including the SPIGOT and two coherent 
de-dispersion systems) will insure that these are only the first of a long series of 
GBT-based results from the double-pulsar. 
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